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1 Introduction 



The gravitomagnetic field is one of tlie most important predictions of general relativity, 
which is believed to have no Newtonian counterpart PI El Ej , that emerges as a conse- 
quence of mass currents, analogous to the generation of magnetic field by electric (charge) 
currents. Working within the framework of general relativity, Lense and Thirring in 1918 
15 predicted that a spinning massive body like Earth would generate a gravitomagnetic 
field, which in turn would cause a precession of the orbits of its satellites. This effect, 
called Lense-Thirring (LT) effect ,is also known as "dragging of inertial frames" or more 
simply "frame-dragging" as Einstein named it. Recently, frame-dragging or LT-effect is 
reported to have been detected [21 El El by Ciufolini et al. [21 El with experimental re- 
sults differing from the general relativistic predictions by about 20% during 1998-2002 
[21 and bylO% in 2004 fS'. Gravity probe-B ( an ongoing space mission of NASA to 
test general relativity using cryogenic gyroscopes in orbit ) was launched in April 2004 
and aims at measurement of Lense-Thirring effect to about 1% [HI error. The important 
fact is that the analysis of the LAGEOS experimental results [21 IS] and the ongoing data 
analysis of Gravity Probe B experiment [HI are based on the assumption that Newton's 
theory of gravitation has no phenomena analogous to magnetism although Newton's law 
of gravitation has a formal counterpart in Coulomb's law of electrostatics. In view of 
the above facts, the possibility of prediction of a Newtonian gravitomagnetic field is of 
immediate practical importance in the context of the gravitomagnetic effects are now 
becoming amenable to experimental observation. By noting the possible existence of a 
gravitomagnetic field first speculated by J. C. Maxwell in 1865 and later pursued by 
Oliver Heaviside in 1893 [8^ , well before the advent of the relativity theories of Einstein, 
here we report the possibility of making a 'natural' prediction of gravitomagnetic field 
within the framework of Newtonian Physics. The full set of Faraday-Maxwell-type field 
equations of gravity describing Newtonian gravitodynamics has been derived here within 
the framework of Newtonian physics in a most natural way. The results represent a sub- 
stantial advance in understanding Newton's theory of gravity in its entirety and also the 
important problem of the gravitomagnetic phenomena in physics. The conclusions seem 
to have immediate and far-reaching implications for the current research on the theory 
and experiments on gravitomagnetic effects in particular and gravitation in general. 



2 History of the Gravitomagnetic field concept 

Recognizing the striking formal analogy between Newton's law of gravitostatics and 
Coulomb's law of electrostatics, J. C. Maxwell [7, , in one of his fundamental works on elec- 
tromagnetism, turned his attention to the possibility of formulating the theory of gravity 
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in a form corresponding to the electromagnetic equations. In this attempt, Maxwell con- 
sidered that the potential energy of a static gravitational configuration is always negative 
but he guessed that this should be re-expressible as an integral over field energy density, 
which being the square of the gravitational field ( by electromagnetic analogy ) is positive. 
Because of this puzzle. Maxwell did not work further on the topic. However, we will see 
a solution to Maxwell's puzzle in this work. 

Nonetheless, in 1893 O. Heaviside pursued Maxwell's attempt further and wrote 
down the full set of Lorentz-Maxwell type equations for gravity, very much analogous to 
the corresponding equations in classical electrodynamics, by virtue of his power of specu- 
lative thought. Heaviside's field equations implied the existence of gravitational waves in 
vacuum, so he considered that the propagation velocity of gravitational waves in vacuum 
might well be the speed of light in vacuum. He also explained the propagation of energy 
in a gravitational field, in terms of gravitoelectromagnetic Poynting vector, even though 
he (as Maxwell did ) considered the nature of gravitational energy is a mystery. Lacking 
experimental evidence of gravitomagnetic effect and for some other reasons , he did not 
work further. As per a report by McDonald (see McDonald in jS]), surprisingly Heaviside 
seemed to be unaware of the long history of measurement of the precession of Mercury's 
orbit. 

The formal analogy was then explored by Einstein [H], in the framework of General 
Relativity. Any theory that combines Newtonian gravity together with Lorentz invariance 
in a consistent way, must include a gravitomagnetic field, which is generated by mass 
current. This is the case, of course, of General Relativity. May be Einstein had not seen 
Heaviside's field equations when he was working on his relativistic theory of gravity. Had 
Einstein seen Heaviside's field equations, his remark on Newton's theory of gravity would 
have been different than what he made before the 1913 congress of natural scientists in 
Vienna pUj.viz.. 

"After the untenability of the theory of action at distance had thus been proved 
in the domain of electrodynamics, confidence in the correctness of Newton's 
action-at-a-distance theory of gravitation was shaken. One had to believe 
that Newton's law of gravity could not embrace the phenomena of gravity in 
their entirety, any more than Coulomb's law of electrostatics embraced the 
theory of electromagnetic processes." 
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3 Derivation of the field equations of Newtonian Grav- 
itodynamics 

In Galileo-Newtonian physics, the source of gravitational field is mass and mass is a 
conserved quantity - any flow of mass must come from some supply, which we know to be 
true. The generation of gravitational field is expressed by Gauss's law of gravitostatics : 

V • 4 = -47rG'p (1) 

— * — * 

where Eg is the gravitational field intensity, p (= mass density) is the source of Eg and 
G is Newton's gravitational constant. The law of conservation of mass is expressed by 
the continuity equation : 

l^-V.J ,2) 

where j is mass current density. Assuming the validity of the two laws (1) and (2) ,we 
now proceed to find the gravitational effects associated with mass currents. We know that 
there exist physical systems where the two laws (1) and (2) are at play simultaneously 
or co-work peacefully. To explore the gravitational behavior of such systems, we take the 
time derivative of (1) to obtain the equation 

1 dEg , dp 



+ ^ = (3) 



inG dt dt 
Now by using (2) in (3), we obtain the equation : 

The quantity inside the parenthesis of (4) is a vector whose divergence is zero. Since 
V ■ (V X y4) = is true for any vector A , the vector inside the parenthesis of (4) can be 
expressed as the curl of some other vector, say Hg . Mathematically speaking, the eq. (4) 
admits of two solutions, viz., 

which is formally analogous to Ampere- Maxwell's Law of classical electrodynamics,viz., 

V X ii" = Je + — (in S.I. units) (6) 
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where H is the magnetic field, je is electric current density and D is the displacement 
vector. In empty space D = cqE , where eo is the permitivity of empty space and E is 
the electric field. One of the solutions (5) would correspond to the reality and the other 
might be a mathematical possibility having no or new physical significance in which we 
are not presently interested in. For our present purpose, we have to choose one of these 
two solutions that may correspond to the reality. Which one to choose ? We can answer 
this question by following the rule of 'study by analogy'. To follow this rule, let us write 
down the Gauss's law of electrostatics in the form : 

V-E = (7) 
Co 

where pe is the charge density, E is electric field intensity,eo is the permitivity of empty 
space. Eq.(7) is analogous to Eq.(l),with only a difference in the sign of the source 
functions. From (6) and (7) we find that pe and je have got the same sign. So by analogy 
we infer that in case of gravity p and j should have the same sign and the signs of je 
and j should be opposite. So by this logic, our most natural choice for the real solution 
of (4) is 

The comparisons of (1) with (7) and (6) with (8) suggest us to introduce or define the 
following terms in gravitation,viz., i^^ as the gravitoelectric (GE) field, eg^ = as 
the gravitoelectric permitivity of empty space, Dg = e^gEg as the gravitational dis- 
placement vector in empty space, Hg as the gravitomagnetic field, and name the Eq.(8) 
as the gravitational Ampere-Maxwell Law of Newtonian gravitoelectromagnetic (GEM) 
theory or Newtonian gravitodynamics. The appearance of the sign difference of the corre- 
sponding source functions in (l)and (7), (6) and (8) imphes the characteristic dissimilarity 
between the two fundamental interactions of nature - which are opposite in nature. For 
instance, in electromagnetism like charges repel and unlike charges attract under static 
conditions, but under dynamic condition the nature of the interaction gets reversed - like 
electric currents (i.e. parallel currents) attract and unlike (i.e. anti-parallel) currents 
repel. In case of gravitation, we have the opposite situation, viz., like masses attract [ and 
if negative mass [IT] exists,the unlike masses should repel ] under static condition, and 
by the nature analogy between gravitational and electrical phenomena, we would have a 
reversed situation in the dynamic condition, viz., like (i.e. parallel) mass currents would 
repel ( as a form of anti-gravity ^2]) and unlike (i.e. anti-parallel) mass currents would 
attract. The deep analogy between electrical and gravitational phenomena in the frame 
work of general relativity has earlier been discussed by R. L. ForwordfT? , R. Wald jTSj, 
V. B. Braginsky,C.M. Caves and K.S. Thorne |14j. One of the implications or predictions 
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of gravitational Ampere-Maxwell Law (8) is that mass currents would generate gravito- 
magnetic field in accordance with what we may call the gravitational Ampere's Law of 
gravitomagnetostatics,viz., 

^ ^ ^ dE 

V X Hg = - j ( when ^ = in (8) ) (9) 
which is the gravitational analogue of Ampere's Law in magnetostatics,viz., 

- - - (9D 

W X H = je (when 7^ = in (6) ) (10) 

To understand the significance and the implications of the 2nd term on the right hand side 
of (8),let us note the significance and the implications of the term ^ (= Cq^) in Ampere- 
Maxwell's Law (6). This term was added to Ampere's Law (10) by Maxwell who called it 
the displacement current. The displacement current is of crucial importance for rapidly 
fluctuating fields. Without it there would be no electromagnetic wave PSl- Analogously 

the gravitational displacement current (= ^og^§f-) is of crucial importance for rapidly 
fluctuating gravitational fields. Without it there would be no gravitational wave. Because 
of these implications of gravitational Ampere-Maxwell Law (8), it is natural to assume the 
existence of gravitational waves in empty space where the law still stands even if j = . 
Thus if gravitational waves exist in empty space, then the electromagnetic analogy we 
are uncovering suggests that the wave equation for the fields Eg and Hg in empty space 
should stand as 

V^-^.-^#-0 (11) 

where Cg is some finite speed of propagation of gravitational waves, which is to be de- 
termined by experiments. To explore this possibility, let us take the curl of (8) to obtain 
the equation 

^'^^ + A^Ft^^""^'^ = Vxj + V-(V-#,) (13) 
In empty space (13) reduces to the wave equation (11), provided the following relations 

AttGOH, , 

^^^« = — IT 

V ■ #3 = (15) 
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V X J = (16) 

hold good in empty space. The relation (16) is certainly true in empty space where j = . 
But the relations (14) and (15) have to be accepted if the existence of gravitational waves 
in empty space is accepted. We now recognize that we have arrived at the full set of 
Faraday-Maxwell-type field equations of gravity which may be represented by 

V ■ Eg = -AnGp = -p/eog, by defining eog = l/47rG (17) 

V X Bg = -fiogj + {l/cg^){dEg/dt), bydefiuiug fiog=AnG/cg^ (18) 

V-Bg = (19) 
V X Eg = -dBg/dt (20) 



where we have defined Bg = {A-nG / Cg^)Hg = fJ^ogHg as the gravitomagnetic induction 
field in empty space and fiog as the gravitomagnetic permeability of empty space, in 
analogy with the electromagnetic case where magnetic induction B in empty space is 
defined by -B = fioH , /io is the magnetic permeability of empty space. This definition 
of /iog ensures the relation Cg = (eo^/iog)"^'^^ in complete analogy with its electromagnetic 
counterpart : c = (eoyUo) . As in electrodynamics, /iog would be the coefficient that 
would determine the strength of gravitomagnetic interaction or effects. Assuming Cg = c, 
a typical value of /ioc, may be estimated at /iog = 9.33 X 10~2^A^.sV%^- This gives us 
a glimpse of the order of magnitude of fiog , for we know not yet the exact value of Cg 
to calculate fiQg . It is this order of smallness of the the value of /log that makes the 
strength of gravitomagnetic effects or interaction very negligible. From (18) it is now 
clear that very large mass currents or very rapidly fluctuating fields are required for 
production of gravitomagnetic field Bg of appreciable strength. Therefore any search for 
gravitomagnetic effects should search for physical systems or processes where the mass 
current density is very large or the gravitational field fluctuation is very rapid. Such 
conditions or situations are available in certain astrophysical systems or events. The 
smallness of the value of fiog also explains why the existence of gravitomagnetic effects 
had escaped the human detection in spite of the long history of its study. By the way, it 
is to be noted that in 1893 Heaviside [H] speculated exactly these equations (17-20) with 
different notations and definitions of the various quantities involved. 
To complete the dynamical picture we ask : What replaces the equation F = mEg to 
describe the force on a particle of mass m , when that particle moves with some velocity v 
in given gravitoelectric and gravitomagnetic fields Eg and Bg ? Because of the rich and 
detailed correspondence between electrical and gravitomagnetic phenomena uncovered 
above, one may at once suggest ( as Heaviside did ) the force on mass m is now 

F = mEg + mvx Bg (21) 
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where v is the velocity of mass m in this expression, in analogy with the Lorentz force 
Law of electrodynamics : 

F = qE + qv X B (22) 

where the symbols have their usual meanings. It is to be noted that the gravitational 
Lorentz force law (21) can also be derived from other Newtonian assumptions chosen by 
Feynman and Dyson ^Hj and following the path of Feynman's proof of Mawwell's equa- 
tions as reported by Dyson with his editorial comment on the derivations. The other 
way to arrive at the above Lorentz-Maxwell-type equations of gravity may be that of the 
Newtonian way Schwinger et al.jT2j have chosen to infer the Lorentz-Maxwell equations 
of classical electrodynamics. From all these variant approaches which lead to the same 
equations as detailed above, we are convinced that the set of four equations (17)-(20) 
would form the basis of all Newtonian gravitoelectromagnetic phenomena . When com- 
bined with the gravitational Lorentz force equation (21) and Newton's 2nd law of motion, 
these equations would provide a complete description of the dynamics of the interacting 
massive particles and gravitomagnetic fields. 

To compare the field equations (17)-(20) with those predicted by general relativity 
(GR), we note the following approximations to the Maxwell-type field equations of GR 
in the parametrized-post-Newtonian (PPN) formalism [Tl] .which as per the present defi- 
nitions and notations may be written as 

V -Eg = -AixGp (23) 

Vx4 = -^ (24) 

(25) 

Vx5,^^-A. + -A.j^-^p.^+^.:^j (26) 

where Ai and A2 are PPN parameters, p is the density of rest masses in the local frame 
of the matter, v is the ordinary (co-ordinate) velocity of the rest mass relative to the 
PPN co-ordinate frame. In general relativity Ai + iAa) = 4 and so Eq.(26) can be 
rewritten as 

- - IQtiG ^ 4 dE„ 

VxB,,^-— p„+-.^ (27) 
In empty space (where p = 0), these field equations reduce to the following equations: 

V ■ 4 = (28) 




7 



Vx4 = -^ (29) 
V-Bg = Q (30) 

Vxi.^l.f (31) 

As per our previous analysis, these Maxwell-type equations of GR imply that 

eog,GR = fXog,GR = 16tt G / , (32) 

which in turn imply that in the low velocity and weak field approximation of GR the 
speed of gravitational wave in empty space is Cg^cR = {^og,GR^Qg,GR) = c/2, a result 
not expected from a Lorentz covariant theory of gravitation. This value of Cg = c/2 can 
also be inferred from the wave equations that follow from the above equations (28)-(31) 
by taking the curl of (29) and utilizing Eqs.(28) and (31) we get the wave equation for 
the field Eg in empty space as 

V^-4-A-^ = (33) 



c, 



g 



where Cg = c/2. Similarly the wave equation for the field Bg can be obtained by taking 
the curl of Eq.(31) and utilizing Eqs.(29) and (30) : 

where again we get Cg = c/2 . However, Peng in |^ has discussed a set of Maxwell-like 
equations that arise in the slow motion, weak field limits of Einstein's field equations 
which in the present notation agrees with Newtonian gravitodynamic equations (17-20) 
with Cg = c [ini 120] • In [201 Peng's equations have been utilized to predict the quantiza- 
tion of planetary revolution speeds which matches with the experimental data. It is to be 
noted that Ciufolini's theoretical analysis of LAGEOS data is based on the gravitomag- 
netic field Bg that can be obtained from Eq.(27) which differs from Peng's corresponding 
equation by a factor of 4 [THl HOI 1201 which is responsible for the above result Cg = c/2 . 



4 Immediate theoretical consequences of Interest 

Newtonian gravitodynamics is very much analogous to Maxwell's electrodynamics as re- 
vealed by the form of its equations. Therefore gravitational phenomena very much anal- 
ogous to those of electromagnetic theory are not surprising to be revealed by this theory. 
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However few concepts and results of unconventional nature and importance may be dis- 
cussed as under. Historians tell us that Newton was quite unhappy over the fact that 
his law of gravitation implies an act ion- at- a- distance interaction over very large distances 
such as that between the Sun and the Earth. But he was unable to resolve this prob- 
lem To resolve Newton's problem within the framework of Newtonian physics and in 
terms of potential functions of Newtonian gravitodynamics, we note that the homogeneous 
equations (19) and (20) admit of the solutions : 

Bg=VxAg, 4 = - V ■ - dlg/dt (35) 

where $g and Ag represents respectively the gravitational scalar and vector potential of 
this theory. These potentials satisfy the inhomogeneous wave equations : 

V' • "^^^ - ^ ■ ^ = 4^^/^ = P/^o, (36) 
^2 7 1 d'^A 47rG'^ 

V^-^. -—2-^ = = Po,J (37) 

^ 9 

if the gravitational Lorenz gauge condition 



is imposed. These will determine the generation of gravitational waves by prescribed mass 
and mass current distributions. Particular solutions (in vacuum) are 



<S.,(f,t)^-G j ^j^^^dv' (39) 

4(.^0^-^/ff^<i.' (40) 

where t' = t — \ r — r' \ / Cg is the retarded time. These are called the retarded potentials. 
Thus we saw that retardation in Newtonian gravity is possible in flat space and time in 
the same procedure as we adopt in electrodynamics. Had Einstein seen these possibilities, 
his confidence in the correctness of Newton's theory would not have been shaken and his 
approach to relativistic gravity would have been different with the implication that the 
history and the character of gravitation might have been different from what we know 
today. Now let us come to Maxwell's puzzle over the gravitational field energy density. 
Actually this puzzle is not a real puzzle, but a guess-work of Maxwell, since the actual 
calculation done by electromagnetic procedure yields the energy density of gravito-electric 
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(i.e. the electric-type component of gravity) and gravitomagnetic (i.e. the magnetic-type 
component ) field , respectively as 

1 - - 1 - - 

(i) Uge = - -tOgEg ■ Eg, {ii) Ugm = - Bg ' Bg (H) 

where e^g = l/AnG and /zog = 4:TTG/cg'^ and the total field energy density is given by 
a sum of the above two,i.e. 

Ufield = Uge + Ugm (42) 

For a particle at rest, i.e., in gravitostatics, the only contribution to its gravitational field 
energy is that due to the gravito-electric field. In gravitostatics ,it is easy to compute 
the gravitational or gravito-electric self energy of a sphere of radius R and mass M with 
uniform mass density by using Eq.(41i), which comes out as : 

U, ^ -1.^,1 E,H.rUr ^ (43) 

The result (43) is in complete agreement with the Newtonian result. It is to be noted that 
Visser [22] used exactly this definition of gravitational field energy density in his classical 
model for the electron. Such a definition of the field energy of gravity has the advantage 
of describing the correct nature of gravitation on quantization because in analogy with 
electromagnetic theory, the present theory will eventually lead to a gauge field of spin 1 
and the spin 1 gauge fields having positive and definite field energy, on quantization, as 
we know lead to a repulsive force field for identical charges of such fields. It is due to this 
reason Gupta [211 Feynman suggested rejection of any spin 1 gauge theory of 
gravity with field energy being positive and definite as such fields do not account for the 
observed nature of gravitational interaction. Since general relativity speaks of spin 2 tensor 
theory of gravity and Newtonian theory as developed here implying spin 1 vector theory 
of gravity very much analogous to the electromagnetic theory, one of the fundamental 
questions raised in the Dicke Framework (see in [21]) : 

" What types of fields, if any, are associated with gravitation - scalar fields, 
vector fields, tensor fields, ?" 

requires immediate attention to understand the riddle of gravitation better. The very 
possibility of a Newtonian gravitoelectromegnetic (GEM) theory also raises another im- 
portant question : " Does general relativity obey the correspondence principle ?" When 
it does, (as per a comparison of Peng's [TH] work with the present one assuming Cg = c ) 
do the predictions of gravitomagnetic phenomena agree with experimental results 01^? 
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5 Analysis of Earth Satellite (LAGEOS) Data 



As an immediate application of Newtonian gravitodynamics, let us consider the analysis 
of Earth satellite (LAGEOS and LAGE0S2) data [21IS]-In the slow motion and weak field 
limit of GR, the gravitomagnetic induction field -Bg.cH associated with the spin-angular 
momentum J of a spherical spinning body ( such as the Earth ) can be estimated from 
(27) following the standard electromagnetic procedure of estimation of magnetic field 
generated by localized current distributions ^^l- The result is that 

Bg^GR = 2G[Jr^ - 3(J - f)-r]/(cV5) (44) 

where f is the position vector of the field point from the center of the Earth. For a 
satellite orbiting the Earth having a gravitomagnetic field presumed at that given by 
(44), the Lense-Thirring (LT) nodal precession is estimated jT] at 

^GR _ 2GJ 
^^^^ " c2a3(l - e2)3/2 

where a and e respectively represents the semi-major axis and the eccentricity of the 
satellite orbit. In 2004 Ciufolini and Pavlis [21 measured the Lense-Thirring nodal pre- 
cession at 

i^experiment ■_^GR 

Q^^ = (0 ■ 99 ± ■ 10)^^^ (46) 

But the above interpretation has not taken care of the contribution arising out of the 
Newtonian gravitomagnetic field, which as we saw here a theoretical possibility. By the 
same procedure, one can estimate the Newtonian gravitomagnetic field of Earth at 

Bg,Neu, = G[Jr^ - 3(J • f)-r\/{2cy) (47) 

and this field would then cause an LT-type nodal precession at 



- 2c,2a3(l-e2)3/2 - [^J (48) 

To see the implications of (48) we consider the following two cases of inter est, viz., 

(a) if Cg = c, then flj^j, = i fi^^ = 25% of ,and 

i^New i^GR i^GR 

(b) if Cg = c/2 , then fi^^ = = 100% of . 

It is to be noted that the exact value of Cg in the slow motion and weak field approxi- 
mation of GR is not unique as seen in sec. 3 of this paper. Because of the uncertainties 
in the exact value of Cg and the possible existence of a Newtonian gravitomagnetic field. 
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the interpretation of the LAGEOS experimental result [S] as a pure general relativistic 
effect with 99% claimed accuracy of GR is doubtful. It may be noted that the 2004 
claim of Ciufolini and Pavlis ^ had been earlier suspected by lorio |23 EH] from other 
possible sources of gravitational error which of course have already been addressed in the 
authors' 2005 work ^ as completely unfounded. Nonetheless, lorio [22] still suspects the 
correctness of the latest results However the most likely Newtonian gravitomagnetic 
error introduced by the result (a) considered above for the case Cg = c seems to be large 
enough to be of some concern let alone the possibility Cg = c/2 leading to the result (b) 
which appears to be more serious one . 



6 Concluding Remarks 

The possibility of a Newtonian gravitomagnetic field can not be ruled out in principle 
and therefore should be taken care of in the current and future measurements of gravit- 
omagnetic phenomena. This possibility adds to our understanding of Newtonian gravity 
in its entirety and raises some fundamental questions of importance in the description 
of gravitational phenomena which may deserve certain attention. The analysis and the 
interpretation of the earth satellite (LAGEOS and LAGE0S2) data [21E1 which ignore(s) 
the Newtonian gravitomagnetic contribution to the observed effects may be reconsidered 
or re-examined for an unambiguous explanation of the experimental observation. 
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